Stickiness in mushroom billiards 
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We investigate dynamical properties of chaotic trajectories in mushroom billiards. These billiards 
present a well-defined simple border between a single regular region and a single chaotic com- 
ponent. We find that the stickiness of chaotic trajectories near the border of the regular region 
occurs through an infinite number of marginally unstable periodic orbits. These orbits have zero 
measure, thus not affecting the ergodicity of the chaotic region. Notwithstanding, they govern 
the main dynamical properties of the system, fn particular, we show that the marginally unstable 
periodic orbits explain the periodicity and the power-law behavior with exponent 7 = 2 observed 
in the distribution of recurrence times. 
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The stickiness of chaotic trajectories in Hamilto- 
nian systems, characterized by long tails in the re- 
currence time statistics, is usually associated with 
the presence of partial barriers to the transport in 
the neighborhood of hierarchies of Kolmogorov- 
Arnold-Moser (KAM) islands. However, as we 
show, these hierarchical structures are not nec- 
essary for the occurrence of stickiness. Here 
we study mushroom billiards which are an- 
alytically solvable systems without hierarchies of 
KAM islands, and we show that these systems 
present stickiness due to the presence of one- 
parameter families of marginally unstable peri- 
odic orbits within the chaotic region. 



I. INTRODUCTION 

An important property of the dynamics in Hamiltonian 
systems with mixed phase space, where regions of regular 
and chaotic motion coexist, is the stickiness (2) of chaotic 
trajectories near the border of KAM islands. These 
chaotic trajectories experience long periods of almost reg- 
ular motion, which strongly influence global properties of 
the system, such as transport (|3|) and decay of correla- 
tions (2). Even though no universal scenario exists 10), 
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stickiness in Hamiltonian systems is in general related to 
topological properties of invariant structures in the phase 
space I n the widely studied case of Hamiltonian 

systems with hierarchies of infinitely many KAM islands, 
there has been an intense debate about the possible uni- 
versality of the scaling exponents in the recurrence time 
statistics (RTS) |1 II II E @). The difficulty in that 
case is partially due to the complexity of the invariant 
structures in the phase space. It is thus of substantial 
interest the study of Hamiltonian systems with simpler 
phase space, where the RTS can be studied analytically. 

Recently, Bunimovich introduced a new family of 
Hamiltonian systems, the so called mushroom billiards, 
which have the remarkable non-generic property of hav- 
ing a phase space with a single KAM island and a single 
ergodic chaotic region. In FIG.^ we show an example of 
mushroom billiard consisting of a semi-circle (hat) placed 
on top of a triangle (foot). In the configuration space, the 
border between the island and the chaotic sea consists of 
the orbits in the semi-circular hat of the billiard that are 
tangent to the circle of radius r in FIG. [IJl The regu- 
lar region (KAM island) is composed of the trajectories 
that do not cross this circle, remaining forever in the 
semi-circular hat. The chaotic component consists of the 
complementary set of trajectories, those that cross the 
dashed circle and may visit the foot of the mushroom. 
Except for zero measure sets, which are discussed in de- 
tail in this article, all trajectories in the chaotic region in- 
deed visit the foot of the mushroom. The stickiness close 
to the border of the chaotic region does not depend on 
the specific format of the foot as far as the foot provides 
a chaotic injection of the trajectories into the hat (where 
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FIG. 1 A mushroom billiard with triangular foot, (a) Con- 
figuration space, (b) Phase-space representation of the semi- 
circular hat, where x is the normalized position and 6 is 
the normalized reflection angle. The parameters are r/R = 
0.6125 and h/R = 1.5, what implies 6* = ±1/6. 



the stickiness occurs). The coordinates used to describe 
the phase space are the normalized position x G [0, 1] 
along the boundary of the billiard and the normalized 
angle 9 6 [—0.5,0.5] with respect to the normal vec- 
tor right after the specular reflection. With respect to 
these coordinates, the border between the regular and 
the chaotic regions for trajectories in the semi-circular 
hat is at 6*(r/R) = ±sin~ (r/-R), as shown in Fig. ^d. 
The important control parameter is the ratio r/R and 
the time is counted as the number of reflections. Precise 
experimental realizations of billiards, e.g., in microwave 
cavities and atom optics l)lOt) . have generated additional 
motivation for the investigation of the dynamics of these 
billiards. 

In this article, we study the properties of stickiness 
in mushroom billiards. Two important features are ob- 
served in the RTS: (i) periodicities in the sequence of 
recurrence times and (it) a power-law distribution with 
scaling exponent 7 = 2 in the limit of long recurrence 
times. Both effects are consequence of the presence of 
one-parameter families of marginally unstable periodic 
orbits (MUPOs) inside the chaotic region, which are 
studied here in detail. While the specific properties of 
the former effect are strongly dependent on the partic- 
ular value of the control parameter r/R, the latter is 
invariant under changes of r/R. 

The article is organized as follows. In Sec. [HI we mo- 
tivate the problem with numerical observations of sticki- 
ness in mushroom billiards. The general characterization 
of MUPOs in mushroom billiards is presented in Sec. IIIII 
In Sec. IIVI MUPOs are used to describe the effects (t) 
and (it) in the stickiness of the chaotic trajectories. A 
summary of the conclusions is presented in the last sec- 



tion. 



II. NUMERICAL OBSERVATION OF STICKINESS 

In this section, numerical observations of the sticki- 
ness in mushroom billiards are reported. We quantify the 
stickiness through the RTS, as follows. The whole foot 
of the billiard is chosen as the recurrence region. A typ- 
ical trajectory is initialized within the recurrence region 
and followed for a long time. The time T the trajectory 
takes, after leaving the recurrence region, to return to it 
for the first time is recorded. Through further iterations 
of the same trajectory, an arbitrarily long sequence of 
recurrence times, distributed according to P(T), can be 
obtained. The RTS is then defined as 

00 N 

T—T 

where TV is the total number of recurrences and N T is the 
number of recurrences with time T > t. For long times, 
in hyperbolic chaotic systems, the RTS decays exponen- 
tially ((111), while in systems with stickiness, it decays 
roughly as a power law (@;|1), 

Q(r) cx r-\ 

where 7 > 1 is the scaling exponent (see also Ref. |l2j) '). 
In our numerics, we approximate Eq. Q with a finite, 
statistically significant number N of recurrences. 

Our main observations about the RTS in the mush- 
room billiard are illustrated in FIG. |21 and can be sum- 
marized as follows: 

(i) The recurrence times T for which recurrences are 
observed appear in a very organized way: times 
without a single recurrence (P(T) = 0) are period- 
ically interrupted by times with a high recurrence 
time probability. The period to between successive 
times with positive probability (P(T) > 0) strongly 
depends on the control parameter r/R and may 
change over large intervals of time T. In particular, 
as shown in the inset of FIG.[21for recurrence times 
in the interval 50 < T < 150, this period is to = 5 
for r/R = 0.5 and to = 11 for r/R = 0.6125. For 
longer recurrence times, higher periods may coexist 
with short periods. 

(it) The overall behavior of the RTS Q(t) shown in 
FIG. [3 presents a clear power-law tail with expo- 
nent 7 = 2, independently of the parameter r/R. 

These observations are fully explained through the 
analysis of MUPOs present in the chaotic component of 
the billiard. The existence and properties of these orbits 
are discussed in the next section. 
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FIG. 2 RTS of the mushroom billiard for various choices 
of the control parameter r/R. The results are consistent 
with a power-law tail with exponent 7 = 2. The distribu- 
tion for r/R = 0.75 was shifted vertically upward by one 
decade for clarity. Inset: the distribution P(T) in the inter- 
val 50 < T < 150 for r/R = 0.5 (to = 5) and r/R = 0.6125 
(to = 11). 

III. MUPOS INSIDE THE CHAOTIC REGION 

The chaotic trajectories in the mushroom billiard al- 
ways visit the foot of the mushroom while the regular 
trajectories are confined to the hat of the billiard. The 
integrability of the latter trajectories is based on the 
conservation of the reflection angle 9 for collisions in 
the semi-circular hat. The novel aspect shown here is 
that there are angle-preserving periodic orbits inside the 
chaotic region that never visit the foot of the mushroom, 
and they form one-parameter families of MUPOs. We use 
the acronym MUPO to refer specifically to marginally 
unstable periodic orbits belonging to this class of orbits. 
They have zero Lyapunov exponent and, in contrast with 
elliptic points, their eigenvalues are real and have mod- 
ulus one. In this section, we study the stability and dis- 
tribution of these MUPOs based on the analysis of an 
equivalent open circular billiard. 

A. Open circular billiard 

A convenient way to visualize the MUPOs is to con- 
sider a circular billiard of radius R, as depicted in 
FIG.|3i, which has the property that trajectories are con- 
sidered to escape when they hit the horizontal straight- 
line segment of length 2r in the center of the billiard 
(hereafter referred to as the hole). The equivalence be- 
tween the two billiards is based on the application of 
the image construction "trick" to the horizontal seg- 
ments in the hat of the mushroom billiard and on the 



independence on the shape of the foot (fjj) . The co- 
ordinates of the circular billiard are the reflection an- 
gle 9 G [—0.5, 0.5] with respect to the normal vector and 
the position of collision in the circumference, given by the 
angle cj) € [— tt, tt], as indicated in FIG. 13^. The time is 
again measured as the number of reflections at the bor- 
der of the billiard. This introduces a minor difference 
between the two billiards since in the mushroom billiard 
one counts the reflections at the horizontal segments of 
the hat. Nevertheless, the dynamics of the trajectories 
in the open circular billiard is equivalent to that of tra- 
jectories in the semi-circular hat of the mushroom bil- 
liard, where stickiness occurs and where the MUPOs are 
located. Geometrically, the MUPOs are the periodic or- 
bits of the open circular billiard that cross the circle of 
radius r but that do not hit the hole in the center of this 
circle. Examples of MUPOs are shown in Figs. |21 and [S] 
for the parameters r/R — 0.5 and r/R = 0.6125. 
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FIG. 3 (a) The open circular billiard for two different lengths 
of the hole: r/R = 0.5 and r/R — 0.6125. Diamonds (♦) cor- 
respond to a periodic orbit (q = 3, r\ = 1) and squares (■) 
to a periodic orbit (q = 4,n = 1). (b) Phase-space rep- 
resentation of (a). Circles (•) correspond to a periodic or- 
bit (q = 7, n = 2), which is studied in detail in FIG. |S] 
and the other symbols are the same as in (a). The four 
horizontal lines represent the border 9* between the chaotic 
and regular regions. The first escape regions for r/R = 0.5 
and r/R — 0.6125 are the areas limited by the dashed and 
dotted curves, respectively. 

We now describe the necessary and sufficient condi- 
tions for the existence of MUPOs inside the chaotic re- 
gion of the mushroom billiard. A periodic orbit in the 
open circular billiard is defined by two integer numbers: 
the period q and the rotation number r\ (i.e., the number 
of laps around the center of the circle), where q/rj is an 
irreducible fraction and q > 2r\. The invariant reflection 
angle 9 of this orbit is given by 9 p (q, rf) — ±(q — 2rj)/2q. 

The first necessary condition for the existence of a 
MUPO defined by q and r\ is that the corresponding tra- 
jectory crosses the circle of radius r, meaning that the 
trajectory would be inside the chaotic component of the 
original mushroom billiard. Accordingly, the reflection 



angle of this orbit must satisfy 

\8p(q>v)\ - 
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\6*(r/R)\. (2) 



The second necessary condition for this orbit to be a 
MUPO is that its trajectory does not hit the hole. In 
order to study this condition in the phase space of the 
open circular billiard, we analyze the points with reflec- 
tion angle 9 that hit the hole in one time step. These 
points define what we call the first escape region, whose 
width wc denote by A r (0) (see FIG.|3Jd). The trajectory 
does not hit the hole if all the q periodic points of the 
orbit are outside the first escape region. With respect to 
the coordinate <p, the distance between two neighboring 
periodic points is constant, namely 2ir/q. The second 
condition can thus be written as 



2tt 



>jA r (0 p ), j 



1 if q is even, 

2 if q is odd. 



(3) 



The factor j = 2 for the odd-period periodic orbits comes 
from the 2n periodicity of the points of these orbits in 
opposition to the it periodicity of the escape region. 

To calculate A r (9), we have to determine the borders 
of the first escape region (FIG. EJd). For a given posi- 
tion 4> in the circumference of radius R, the angles 9 of 
the trajectories that first hit the hole are limited by the 
angles 6' ± given by 
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The width of the first escape region is then given by 

A r (0 p ) =cf>t- 4 , (5) 



where 



> 02 are the two solutions of the equa- 
tion ± = 9 p . Observe that, because of the symme- 
try 9 — » —9, we have §\ — 4>\ = — ^ . 

Gathering all these, we have that the conditions ex- 
pressed in the Eqs. J2J and @ are not only necessary 
but also sufficient for an orbit with period q and rota- 
tion number rj to be a MUPO inside the chaotic region 
of the original mushroom billiard. These conditions can 
be translated in terms of the interval of the control pa- 
rameter r/R for which a given periodic orbit (q,ri) is a 
MUPO: 



sin[7r6» p (q, rjj\ < r/R < 



sin [nO p (q, Tj)] 



(G) 



COS [7r/0*£?)] 

with j as in Eq. In Fig. 0] we show the parameters 
for which the orbits up to period q = 20 are MUPOs. 
An efficient procedure to find higher order MUPOs for 
a given parameter r/R is to take r\/q as the convergents 
of the continuous fraction expansion of ^ cos -1 (r/R) and 
verify if they fulfill condition © or, equivalently, Eqs. |j5J 
and J3J. Through this procedure we have found the fol- 
lowing MUPOs for r/R = 0.625: (q,rj) = (7, 2), (698, 
199), (1161, 331), (18341, 5229), (2136146, 609013), and 
(8526243, 2430823). 




FIG. 4 (Color) Intervals of the control parameter r/R for 
which orbits (q, r/) are MUPOs. All orbits with q < 20 are 
shown. 



B. Stability of periodic orbits 

In this section we show that the periodic orbits con- 
sidered in Sec. IIII. Al are marginally unstable, i.e., are un- 
stable but have zero Lyapunov exponent, and that they 
form one-parameter families in the phase space. Con- 
sider small perturbations to one of these orbits. If the 
initial position <fi is perturbed, neglecting the presence of 
the hole, an equivalent orbit rotated in the configuration 
space is always obtained (FIG. |5j). When we consider 
the presence of the hole, some of these orbits escape and 
create gaps in the possible <fi values of the periodic or- 
bits. However, when the (strict) inequalities in Eqs. (0) 
and are satisfied, there remain (open) intervals of <f> 
values in between these gaps. These remaining intervals 
define a one-parameter family of MUPOs with constant 9. 
A small perturbation in the angle 9' = 9 — e, on the other 
hand, leads to an increase <j)' — (j) = 2e per time step in 
the distance from the original unperturbed orbit. This 
shows that the separation increases linearly in time and 
hence that these periodic orbits are indeed marginally 
unstable. Since the marginal stability does not depend 
on the discretization of time, the statistical properties of 
the recurrence time discussed in this article remain valid 
for continuous time dynamics. 

Typically, a MUPO perturbed with respect to (0, <f>) 
is a quasi-periodic orbit in the circular billiard, which, 
without the hole, would pass arbitrarily close to all values 
of (p. Because of this ergodicity with respect to the co- 
ordinate <j), the corresponding perturbed trajectory even- 
tually enters the first escape region and hits the hole, as 
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FIG. 5 Detailed analysis of orbits (q = 7, 77 = 2) for 
r/R — 0.6125. (a) Configuration space, where two orbits 
(q = 7, 77 = 2) that cross the circle with radius r are shown. 
The orbit represented by circles (•) does not hit the hole, 
while the orbit represented by triangles (▼) hits the hole on 
the right-hand side (see amplification), (b) Phase-space rep- 
resentation of the orbits in (a), where it is shown that they 
are respectively outside and inside the first escape region. A 
small perturbation in the reflection angle 9 of the first orbit 
leads to a continuous rotation of the orbit in (a) [horizontal 
drift in (b)] until the trajectory hits the hole [enters the first 
escape region in (b)] . (c) Time evolution of the distance from 
the regular region to a chaotic trajectory that approaches the 
family of MUPOs (q — 7, 77 = 2) in the original mushroom 
billiard. Events with recurrence time T — 463 and T = 17 
are highlighted. The events with large recurrence times are 
associated with approaches to the MUPOs. 



illustrated in FIG.|SJ The distance from the unperturbed 
MUPO grows linearly in time. However, the distance e 
from the corresponding family of MUPOs remain con- 
stant until the perturbed trajectory hits the hole and 
escapes (FIG. 03:). 



C. Distribution of families of MUPOs 

We now use the information about the first escape re- 
gion, Eq. to study how the families of MUPOs are 
distributed inside the chaotic region of the mushroom 
billiard. The MUPOs are the only orbits in the mush- 
room's hat that cross the circle of radius r but never 



visit the foot of the mushroom (Fig. . Equivalcntly, 
in the open circular billiard, these are the orbits that 
cut the circle of radius r but never escape (Fig. I n 
FIG. H3 we show the phase space of the open circular 
billiard for trajectories |6*| < \6*\, associated with the 
chaotic region of the mushroom billiard, and the iVth es- 
cape region, defined by the (N — l)th pre-image of the 
first escape region shown in FIG. |3t>- We note that the 
larger the period of the orbit the closer it may be to 
the border 8*. For example, for the control parame- 
ter r/R = 0.6125 considered in FIG. EJ the border is 
at 9* = 0.2149010... and the orbits (q = 7,77 = 2), 
(q = 235, 77 = 67) and (q = 698, 77 = 199) highlighted 
in the figure are at 9 p = 0.2143857..., 6 p = 0.2148936... 
and 6 P = 0.2148997..., respectively. The widths of the 
escape regions go to zero when the border is approached. 
In Fig. |H| the MUPOs correspond to the points that do 
not belong to any of the nth escape regions, for all n < N 
in the limit of N — > 00. That is, taking the limit N — > 00 
in FIG.H3 all the points outside the escape regions belong 
to MUPOs. Note that a complex distribution of families 
of MUPOs may exist near the border of the island. 
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FIG. 6 (Color online) Escape regions in the the phase space 
of the open circular billiard for r/R = 0.6125. In each panel, 
the iVth escape region is shown in black and the nth escape 
regions for all n < N are shown in gray. In the first row, we 
show the cases N — 2, iV = 3, and N = 8, while in the second 
row we show the case N — 20 and successive amplifications 
for this case. Different symbols correspond to the different 
orbits (q = 7,rj = 2), (q = 235, 77 = 67), and (q = 698, 77 = 
199), from top to bottom, respectively. The figures at the 
bottom right show that only the first and the last of these 
orbits are MUPOs. 
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IV. MUPOS AND STICKINESS 

In this section, we show how the MUPOs character- 
ized in Sec. IIIII affect the stickiness of chaotic trajecto- 
ries in the mushroom billiard and explain the numerical 
results reported in Sec. [H] Qualitatively, the stickiness 
of chaotic trajectories in mushroom billiards can be de- 
scribed as follows. A chaotic trajectory bounces in the 
foot of the mushroom until it enters the semi-circular hat. 
Depending on the reflection angle 9, which is conserved 
for collisions in the hat, the trajectory may be close to a 
MUPO. In this case, it remains close to the family of this 
MUPO a long time before falling again into the foot of the 
billiard (see FIG. Et)- In the phase space, these events 
with large recurrence time are observed as a continuous 
rotation of the trajectory around the regular island, at 
a fixed distance from the island, until it enters the first 
escape region. 

We now show how the MUPOs explain quantitatively 
the two effects observed in the RTS and anticipated in 
Sec. HQ 

(i) Regularity in the distribution P(T). When a 
chaotic trajectory visits the hat of the billiard, the ap- 
proach to as well as the escape from the neighborhood 
of a family of MUPOs takes place in a single time step, 
as illustrated in FIG.OJ by a sharp transition before and 
after the interval of minimum distance to the island. Due 
to the injection of trajectories from the foot of the mush- 
room billiard, the approach to the MUPOs (that requires 
an angle 9 close to 9 P ) can only occur close to some of 
the q periodic points of the MUPO. Actually, it happens 
always near the same position in the mushroom's hat, 
right above the intersection point between the boundary 
of the foot and the bottom part of the hat (i.e., one of the 
two horizontal lines). Similarly, the escape occurs always 
when the trajectory is close to the point of the MUPO 
that lies in the bottom part of the mushroom's hat. Due 
to these constraints, the intervals of time a chaotic tra- 
jectory that approaches a family of MUPOs spends away 
from the mushroom's foot form a sequence that can be 
written as T{ = a — 1 + (q + 2r])i, where i £ IN and a is the 
time between the first collision in the semi-circular hat 
and the collision in the bottom part of the hat (close to 
the hole). The period to between successive recurrence 
times with positive probability, reported in Sec.^ is thus 
related to the period q and the rotation number 77 of the 
MUPO inside the chaotic region according to 

to = q + 2 V . (7) 

The factor 2rj comes from the fact that, in contrast with 
the open circular billiard, in the mushroom billiard one 
counts the collisions in the bottom of the hat. In FIG.0 
we show that for r/R = 0.5 the orbit (q = 3,ry = 1) 
is exactly in the border between the regular island and 
the chaotic region. Through Eq. J7J we obtain to = 5, 
explaining the numerical observation. Analogously for 
the parameter r/R = 0.6125, in FIG. [S] we see the or- 
bit (q = 7, 77 = 2) that implies to = 11, exactly as ob- 



served numerically. Higher-order periodicities are associ- 
ated with the existence of additional families of MUPOs 
in the chaotic shown in Sec. Illl.GI 

(ii) 7 = 2 in the tail of the RTS. In the stability con- 
siderations described in Sec. I11I.BI we noticed that when 
we perturb a MUPO, the distance from the unperturbed 
orbit increases linearly along <j> and remains constant with 
respect to the reflection angle 9. These properties of the 
MUPOs in mushroom billiards are equivalent to those 
observed for MUPOs present in the much simpler case of 
billiards with parallel walls. In the latter case, the MU- 
POs are generated by trajectories bouncing perpendicu- 
larly between parallel walls. Billiards with parallel walls, 
such as the Sinai billiard and stadium billiard, have been 
studied extensively |l3t IhH : Hit H^) . It has been shown 
that the existence of a one-parameter family of MUPOs 
due to parallel walls leads to an exponent 7 = 2 in the 
power-law decay of the RTS {3 L3 0)- Each family 
of MUPOs in the mushroom billiards is equivalent to a 
family of MUPOs generated by the presence of parallel 
walls. The same is true for the family of periodic or 
quasi-periodic marginally unstable orbits always present 
at the border between the regular and chaotic region of 
the mushroom billiard. In the mushroom billiard, more 
than one family of marginally unstable orbits may exist. 
Since each of them leads to the same asymptotic expo- 
nent 7 = 2, irrespective of the control parameter r/R, 
the overall RTS will also have scaling exponent 7 = 2. A 
similar argument can be used to show that the same ex- 
ponent governs the RTS of honey mushrooms ( 1), which 
are systems composed by a finite number of mushroom 
billiards. 



V. CONCLUSIONS 

Mushroom billiards (Q) provide an example of non- 
trivial Hamiltonian systems with a relatively simple 
phase space. In contrast with previously considered 
Hamiltonian systems, which have a hierarchy of infinitely 
many regular islands, the mushroom billiards considered 
here have a single island. We have shown that, in spite of 
the simplicity of the regular region, the chaotic compo- 
nent of the mushroom billiards contains a complex distri- 
bution of marginally unstable periodic orbits (MUPOs). 
These orbits have zero measure and do not affect the er- 
godicity of the system. Nevertheless, the MUPOs are cru- 
cial for the understanding of important dynamical prop- 
erties of the system, as shown here through the analysis 
of the recurrence time statistics. In particular, we have 
shown that these MUPOs lead to an exponent 7 = 2 for 
the asymptotic power-law decay of the recurrence time 
statistics. This result leads to the interesting question 
of whether this scaling exponent is universal within the 
class of Hamiltonian systems with finite number of regu- 
lar islands l)l7l) which, in turn, may provide new insight 
into the problem of universality in Hamiltonian system 
in general. 
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